gusucode.com > MATLAB与控制系统仿真实践 源代码 > MATLAB与控制系统仿真实践 源代码/第十四章/eg14_13.m

    wc2=1.5;
num=20;
den=conv([1 1 0],[1 2]);
G=tf(num,den);
[mag,phase,wcg,wcp]=margin(G);
margin(G)
t1=1/(0.1*wcg);
beta=10;
Gc_lag=tf([t1,1],[beta*t1,1])
G1=G*Gc_lag;
[mag,phase,w]=bode(G1);
mag1=spline(w,mag,wc2);
L=20*log10(mag1);
alfa=10^(-L/10);
t2=1/wc2/sqrt(alfa);
Gc_lead=tf([alfa*t2,1],[t2,1]);
G0=G*Gc_lead*Gc_lag;
figure(2)
margin(G0)
figure(3)
step(feedback(G0,1))